AD-AU2  17*»  AERONAUTICAL  RESEARCH  LABS  MELBOURNE  (AUSTRALIA)  F/6  *0/11 

THE  SIMULATION  OF  LOCAL  STRESS-STRAIN  BEHAVIOUR  WITH  APPLICATIO— ETC f U1 
JUL  79  R  C  FRASER  1 

UNCLASSIFIED  ARL/STRUC-305  ^ 


one  auuiQEf 


ARL- STRUC-REPORT- 38 5 


AR-002-247 


*7  DEPARTMENT  OF  DEFENCE 

^  DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
rH  AERONAUTICAL  RESEARCH  LABORATORIES 

MELBOURNE.  VICTORIA 


STRUCTURES  REPORT  385 


THE  SIMULATION  OF  LOCAL  STRESS-STRAIN 
BEHAVIOUR  WITH  APPLICATION  TO  FATIGUE 


R.  C.  FRASER 


Approved  for  Public  Release. 


©  COMMONWEALTH  OF  AUSTRALIA  1981 

CORY  No  (i  JULY,  I?79 


*  1 

v.*  X 


AR-002-247 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  LABORATORIES 


STRUCTURES  REPORT  385 


THE  SIMULATION  OF  LOCAL  STRESS-STRAIN 
BEHAVIOUR  WITH  APPLICATION  TO  FATIGUE 

ANALYSIS 


by 

R.  C.  FRASER 


o. 


\C 


& 


SUMMARY 

Substantial  improvement  in  fatigue  life  prediction  has  been  effected  by  the  introduction 
of  cyclic  stress-strain  concepts.  A  procedure  that  incorporates  these  in  a  simulation  of 
local  stress  and  strain  from  nominal  load  using  Neuber  notch  analysis  is  described  in 
relation  to  its  application  to  the  fatigue  analysis  of  in-service  toad  histories. 

In  the  present  simulation  the  non-linear  stress-strain  and  Neuber  equations  are 
solved  numerically  rather  than  by  piecewise  approximations.  For  the  long  load  histories 
occurring  in  practice  this  improves  numerical  stability. 
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ABSTRACT 


Substantial  improvement  in  fatigue  life  prediction  has  been  effected  by  the  intro¬ 
duction  of  cyclic  stress-strain  concepts.  A  procedure  that  incorporates  these  in  a 
simulation  of  local  stress  and  strain  from  nominal  load  using  Neuber  notch  analysis 
is  described  in  relation  to  its  application  to  the  fatigue  analysis  of  in-service  histories. 

In  the  present  simulation  the  non-linear  stress-strain  and  Neuber  equations  are  solved 
numerically  rather  than  by  piecewise  approximations.  For  the  long  load  histories  occur¬ 
ring  in  practice  this  improves  numerical  stability.  ^  _ - 
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DISTRIBUTION 


I.  INTRODUCTION 


Expressed  simplislieally,  the  economic  design  of  a  modern  metallic  structure  often  dictates 
that  it  operates  with  a  minimum  amount  of  material  for  a  substantial  life.  The  former  require¬ 
ment  is  usually  determined  by  such  factors  as  capital  cost,  size  or  weight  constraints,  while  the 
latter  is  determined  by  the  need  to  maximise  the  gains  for  which  the  structure  was  designed. 
However,  there  exists  a  fundamental  conflict  between  these  two  design  objectives  in  that  the 
former  inevitably  raises  peak  stresses  and  strains  while  forcing  down  the  life  of  the  structure. 
The  trade-off  between  the  two  has  received  much  attention  in  recent  years  not  only  in  relation  to 
design  but  also  in  response  to  a  need  to  re-evaluate  the  condition  of  old  structures  that  are  being 
called  upon  to  last  longer  than  originally  anticipated. 

In  the  literature  various  methods  have  been  developed  recently  to  enable  one  to  account 
for  the  effects  of  peak  stresses  and  strains  on  fatigue  life.  Of  these  a  technique  originally  proposed 
by  Neuber1  for  static  loading  but  developed  by  Topper  ei  uL-  has  obtained  unique  importance 
for  fatigue.  Based  on  the  notch  analysis  of  shear  strained  prisms  the  method  allows  the  deter¬ 
mination  of  stresses  and  strains  local  to  a  specific  stress  raiser  using  information  on  the  nominal 
stress-strain  state  of  the  member  in  question  in  conjunction  with  data  detailing  the  response  of 
the  metal  to  cyclic  plasticity. 

Two  distinct  approaches  to  the  analysis  of  fatigue  life  using  Neuber  analysis  have  emerged. 
Topper  ei  a!.'1.  Morrow  ct  al .?  and  Wetzel4  have  analysed  successfully  the  fatigue  behaviour  of 
notched  structures  by  using  Neuber  analysis  to  control  the  loading  of  appropriate  smooth 
specimens.  Other  researchers5  however,  have  studied  the  fatigue  behaviour  of  various  struc¬ 
tures  by  usine  Neuber  analysis  to  simulate  the  stress  strain  state  local  to  critical  stress  raisers. 
The  fundamental  difference  between  the  two  approaches  is  that  whereas  the  latter  must  use 
empirical  relations  describing  the  stress-strain  and  damage  response  of  the  material  to  cyclic 
load,  the  former  utilises  an  “exact”  representation  of  this  information  the  specimen. 

However,  the  adv  antages  in  terms  of  cost,  time  and  flexibility  of  a  computer  simulation  over 
a  laboratory  simulation  of  local  stress-strain  cum  fatigue  behaviour  is  such  that  it  is  likely  to 
be  used  increasingly  in  the  future  for  the  life  analysis  of  both  new  and  old  structures. 

The  present  paper  describes  the  initial  development  of  a  computer  program  written  in 
PASCAL  that  is  designed  to  predict  the  life  of  structures  under  in-service  load  histories.  The 
advantage  of  this  program  is  that  the  non-linear  model  equations  arc  solved  directly  without 
further  approximation  of  the  stress -strain  equations.  Over  a  long  sequence  this  improves 
numerical  stability  for  a  surprisingly  small  computing  cost. 


2.  CYCLIC  STRESS-STRAIN  CONCEPTS 

As  a  background  this  section  and  the  next  summarise  the  salient,  well  known  features  of 
local  plasticity  in  fatigue. 10,u'12 

The  inelastic  response  of  metals  to  cyclic  load  can  be  complex  as  demonstrated  in  Figure  1 
and  involves  the  simultaneous  interaction  of  stress  strain  hysteresis  or  memory,  cyclic  softening 
or  hardening  and  the  relocation  of  mean  stress  through  cyclic  creep.  Of  these,  hysteresis  is  the 
most  obvious  and  important  and  can  be  defined  simply  as  a  lack  of  reversibility  when  a  given 
stress-strain  state  or  path  moves  towards  a  tension  or  compression.  Figure  2  depicts  the  response 
of  a  metal  at  a  stress  raiser  to  a  completely  reversed  loading,  and  it  can  be  seen  that  the  local 
stress-strain  state  moves  along  two  distinct  paths  when  cycling  to  or  from  either  the  compres¬ 
sive  or  tensile  situations,  i.e.  the  local  stress  strain  path  does  not  retrace  its  movement  when 
cycling  to  or  from  a  given  state.  The  initial  hystcrisis  loop  that  is  evident  in  Figure  2  does  not 
represent  a  stable  condition  in  most  practical  instances  and  will  alter  as  shown  in  Figure  3  upon 
application  of  further  load  cycles;  when  the  strain  change  per  stress  change  for  each  reversal 
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FIG.  1  REPRESENTATION  OF  NOMINAL  AND  LOCAL  STRESS  -  STRAIN 
BEHAVIOUR  AT  A  STRESS  RAISER  DURING  A  COMPLEX 
LOADING  SEQUENCE 


is  decreasing  locally  then  the  material  is  said  to  cyclically  harden  and  vice-versa  for  cyclical 
softening.  For  most  materials,  after  many  cycles,  the  hysteresis  loop  will  eventually  stabilize, 
and  by  plotting  the  locus  of  tips  of  stable  hysteresis  loops  from  companion  tests  conducted  at 
different  amplitudes  of  local  reversed  stress  (or  strain)  one  obtains  a  steady  state  stress-strain 
curve  of  amplitude.  This  is  analogous  to  and  can  be  compared  with  the  usual  monotonic  stress- 
strain  curve — see  Figure  4.  Cyclically  induced  changes  are  immediately  evident  by  such  com¬ 
parisons  as  shown  in  Figure  5.  A  metal  that  cyclically  hardens  will  have  a  cyclic  stress-strain 
curve  (as  the  locus  of  stable  hysteresis  loop  tips  is  called)  above  the  monotonic  curve  and  con¬ 
versely  if  the  cyclic  curve  is  below  the  monotonic  curve  the  metal  cyclically  softens. 

Hardening  and  softening  are  usually  transient  phenomena  that  occur  early  in  the  fatigue 
life  of  a  metal,  often  quite  rapidly*  (see  Section  2.2).  This  frequently  permits  simulation  of  fatigue 
using  only  the  steady  state  properties  without  affecting  life  predictions  significantly.  The  deter¬ 
mination  of  cyclic  stress-strain  curves  using  fast  incremental  step  strain  tests  is  described  in 


*  Typically,  softening  occupies  30%  of  total  life  while  hardening  is  complete  after  l°„. 
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Reference  11.  Cyclic  stress-strain  data  for  a  variety  of  metals  is  detailed  in  Reference  12.  The 
law  most  commonly  used  to  fit  the  cyclic  stress-strain  curve  is  similar  to  that  used  for  mono¬ 
tonic  curves,  namely 


‘t  !  fp 


where  «t  =  total  strain, 
ee  =  elastic  strain, 

€p  =  plastic  strain, 
a  —■  total  stress, 

E  Young's  modulus, 

A"  -  cyclic  strength  coefficient. 
ri  -  cyclic  strain  hardening  exponent. 

The  value  of  the  cyclic  strain  hardening  exponent  is  around  0T5  for  most  metals  (a  mono¬ 
tonic  strain  hardening  exponent  n  ■  n'  will  indicate  that  the  metal  cyclically  hardens  or  vice 
versa). 

Masing111  has  shown  that  hysteresis  loop  shape  can  be  related  to  the  cyclic  stress-strain  curve 
by  simply  doubling  both  scales.  Figure  6  illustrates  this  geometric  similarity  between  a  local 
stress-strain  path  and  some  segment  from  the  start  of  the  cyclic  stress-strain  curve. 

Thus  the  total  change  of  local  strain  in  terms  of  local  stress-increments  is 


which  relates  the  hysteresis  loop  to  the  cyclic  stress-strain  curve  using  the  co-ordinate  system  of 
Figure  7a,  and  parameters  of  the  cyclic  curve.  Figure  7  b  shows  the  two  curves  joined  as  part  of 
a  loop.  The  starred  quantities  distinguish  the  two  axis  systems. 


FIG.  6  MASING'S  RULE  -  RELATION  BETWEEN  CYCLIC  STRESS-STRAIN  CURVE 
AND  HYSTERESIS  LOOP 
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FIG.  9  COMPLEX  MEMORY  BEHAVIOUR 


2.1  Memory 

The  next  aspect  of  cyclic  plasticity  that  needs  to  be  considered  is  memory.  Figure  8  demon¬ 
strates  this  for  a  metallic  test  specimen  which  has  been  strained  in  tension  from  condition  I 
to  condition  2  and  then  partially  unloaded  and  loaded  again.  In  deforming  from  condition 
3  to  4  the  metal  remembers  the  prior  deformation  and  deforms  according  to  the  initial  response 
along  1  to  2.  Under  the  asymptotic  steady  state  the  shape  of  the  hysteresis  (or  memory)  loop 
between  2  to  3  in  Figure  8  is  still  described  by  (2)  (there  used  twice).  Even  in  compounded  situ¬ 
ations  such  as  in  Figure  9.  where  the  local  stress-strain  state  is  moving  from  point  6  to  7  in 
response  to  the  nominal  load,  change  6  to  7  remembers  three  prior  deformation  conditions. 

2.2  Cyclic  Relaxation  and  Creep 

The  last  phenomena  to  be  considered,  cyclic  relaxation  of  mean  stress  and  cyclic  creep  under 
asymmetric  loading  can  be  demonstrated  by  strain-controlled  testing  in  the  first  instance  and  by 
stress-controlled  testing  in  the  second.  In  Figure  ICu  the  relaxation  of  mean  stress  towards  zero 
in  response  to  cycling  between  fixed  local  strain  limits  is  clearly  evident,  likewise  when  the  stress 
is  controlled  as  in  Figure  10b  the  local  strain  "creeps".  Cyclic  relaxation  and  cyclic  creep 
represent  essentially  the  same  response  (related  to  hardening  and  softening)  under  different 
control  conditions  and  can  be  simulated  by  changing  K'  and  n'  in  (2)  during  cycling.  Material 
test  data  that  permits  the  determination  of  the  forcing  functions  A"(7V)  and  n’(N).  say  to  simulate 
cyclic  relaxation  and  creep  is  sparse  but  it  often  leads  to  no  significant  change  in  life 
predictions.4  r'. 


3.  NEUBER  NOTCH  ANALYSIS 

Neuber's  analysis  of  stress  and  strain  concentration  in  notched  prismatic  bodies  in  anti¬ 
plane  shear  has  resulted  in  a  state  equation  that  is  readily  adapted  to  simulation  of  local  stress- 
strains.  In  his  paper  Neubcr  sb  ws  that  for  this  two-dimensional  situation  the  geometrical  mean 
of  the  stress  and  strain  concentration  factors  under  any  arbitrary  non-linear  stress- strain  law  is 
equal  to  the  Hookean  stress  concentration  factor,  i.e. 
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FIG.  10  REPRESENTATION  OF  CYCLIC  RELAXATION  AND  CYCLIC  CREEP 


At  (A.  AJ»  * 

Topper  et  air  have  transcribed  this  equation  to  relate  elastic  change  in  nominal  stress  to 
change  in  local  stress  and  strain  as  follows: 

A,  (A„  A,)1  2 


or 


AaA<  (A,  ISr-jE. 


(4) 


8 


In  this,  a  and  f  refer  to  local  conditions  while  S  and  e  S  £  are  nominal.  Thus  for  a  given 
change  in  nominal  stress  the  changes  in  the  local  stress  and  strain  will  be  related  by  the  rect¬ 
angular  hyperbola  (4).  Several  researchers1-1-'  ’ '-14  realising  that  small  sharp  notches  have  less 
effect  on  fatigue  life  than  indicated  b>  Kt,  have  suggested  that  in  such  circumstances  a  modified 
factor  Kt  be  used  in  lieu  of  K\  in  (3)  and  (4).  Called  a  fatigue  strength  reduction  factor,  it  accounts 
for  size  effects  and  is  used  as  a  constant  for  a  given  material  and  geometry.  Determination  of  Kt 
is  discussed  in  References  2  and  28  but  here  wc  retain  the  notation  K ,  with  this  option  understood. 

Since  most  service  load  environments  are  recorded  as  nominal  strain  histories  another 
appropriate  form  of  (3)  is 

A<t  At  t.  At-(Ac)-.  (5) 

Equations  (4)  and  (5)  represent  a  stress  strain  relation  determined  mainly  by  geometry.  We  have 
already  mentioned  in  (2)  another  relation  determined  mainly  by  the  material.  Satisfaction  of 
both  these  relations  determines  changes  in  stress  and  strain.  Thus  for  a  nominal  stress  or  strain 
movement  as  in  Figure  11.  the  resulting  local  stress  strain  slate  will  lie  on  the  rectangular 
hyperbola  defined  by  (4)  or  (5).  The  path  along  which  the  local  state  moves  to  that  point  will  be 
appropriate  to  the  characteristics  of  the  material  and  its  prior  loading.  Figure  II  illustrates 
this  use  of  the  Neuber  rule  in  determining  local  stress  strain  movement.  For  the  particular  case 
shown  the  material  response  curve  is  defined  by  the  hysteresis  function.  If  the  material  exhibits 
more  complex  behaviour  such  as  memory  effect  as  shown  in  Figure  12  then  the  local  (and 
overall)  stress-strain  movement  will  be  the  intersection  of  the  Neuber  Rule  and  the  material 
response  curve  that  reflects  the  given  conditions. 


FIG.  12  COMPLEX  LOCAL  STRESS-STRAIN  MOVEMENT  USING  THE  NEUBER  RULE 
(A  SINGLE  MEMORY  INFLUENCE  IS  ILLUSTRATED) 


A  mention  of  the  limitations  of  Equations  (.1)  to  (5)  is  in  order  as  there  are  constraints  on 
their  use.  Firstly  since  (4)  and  (5)  were  derived  assuming  nominal  stress  to  be  always  elastically 
related  to  nominal  strain,  they  cannot  be  used  when  this  is  not  so. 

Secondly,  as  Neuber's  analysis  was  based  on  the  mathematics  of  two-dimensional  antiplane 
shear,  combined  stress  at  a  notch  will  not  be  adequately  described  by  (3)  nor  subsequently 
Equations  (4)  and  (5). 

Finally,  the  equations  apply  only  tu  the  crack  initiation  phase  of  fatigue  life;  where  crack 
propagation  occupies  a  significant  portion  of  the  fatigue  life,  it  must  he  predicted  by  fracture 
mechanics. 

4.  SIMULATION  OF  LOCAL  STRESS-STRAIN 

There  are  several  different  approaches  to  simulating  the  local  stress  strain  response  to  a 
nominal  load  sequence. -tf.ss  Martin  ct  al.  have  developed  a  computer  program  that 
simulates  local  stress  strain  behaviour  using  a  theoretical  model  built  from  combinations  of 
springs  and  Newtonian  frictional  sliders.  Wetzel*  and  Conle1*  later  extended  this  approach 
using  "push  down  lists"  to  control  or  more  specifically  activate  the  use  of  cyclic  stress-strain 
data.  For  example  the  power  law  for  monotonic  or  cyclic  stress  strain  curves  (i.e.  Equations 
(I)  and  (2))  is  approximated  by  a  sequence  of  straight  line  segments  which  are  used  in  construct¬ 
ing  local  stress-strain  movements  according  to  a  vector  of  availability  coefficients  (the  avail¬ 
ability  vector  effectively  stores  the  deformation  history  of  the  material).  Cycle  dependent 
hardening  or  softening  and  relaxation  of  mean  stress  or  creep  can  be  considered  as  small  cycle- 
by-cycle  changes  of  parameters  in  the  cyclic  stress  strain  curve  quotation.  As  functions  K'(n) 
and  n'(.V)  these  can  be  easily  incorporated  into  a  segmented  representation  of  the  stress-strain 
equation  by  altering  the  slopes  of  these  segments  appropriately  .  This  approach  has  been  utilised 
successfully  by  Wetzel  and  Conle  and  is  amenable  to  various  control  conditions,  e.g.  simulating 
local  stress  strain  behaviour  of  specimens  subjected  to  controlled  strain-cycling,  controlled 
stress-cycling  or  Neuber  rule  controlled  cycling  (i.e.  random  loading). 

4.1  Present  Method 

The  simulation  technique  to  follow  differs  from  those  of  the  references  above  as  a  result 
of  a  need  to  process  random,  nominal  s.rain  histories  of  unknown  lengths.  Thus  the  procedure 
had  to  be  fast,  use  Neuber  notch  analysis  as  the  control  condition  and  be  numerically  stable — 
since  Neuber  control  uses  the  existing  local  stress  strain  state  to  determine  the  next  any  inherent 
error  in  such  a  step  is  cumulative.  Consequently  the  determination  of  a  local  stress-strain  state 
should  be  done  as  quickly  as  possible,  preferably  in  one  step;  hence  the  attraction  of  deter¬ 
mining  this  state  directly  from  the  interaction  of  the  appropriate  hysteresis  curve  with  the 
Neuber  Control  hyperbola  as  the  fundamental  operation  of  the  model.  Mathematically  the  state 
is  the  solution  of  (2)  and  (5).  with  appropriate  origins,  for  (A«.  At). 

4.1.1  Relaxation 

Initially  the  simulation  incorporated  cyclic  hardening  softening  and  cyclic  relaxation  creep 
functions  which  modified  A '  and  n  as  appropriate  though  subsequently  these  phenomenon  were 
deleted  from  the  model  as  the  marginal  changes  in  life  predictions  that  accompanied  their  use 
did  not  justify  the  penalties  incurred  w:tri  respect  to  program  speed.  (The  sparsity  of  accurate 
data  that  permitted  quantification  of  the  parameters  was  also  a  problem.)  Comprehensive 
validation  of  simulation  models  by  Wetzel1'  and  Conle  and  Topper"  have  suggested  that  it  may 
unnecessary  to  include  relaxation  in  most  instances. 

4.1.2  Modelling 

Flence  the  basic  simulation  required  only  hysteresis  and  memory  capability. 

For  this  situation  there  exist  only  four  basic  types  of  movement  in  the  local  stress  strain 
state  and  these  are  shown  in  Figure  If  movement  along  the  cyclic  curve  0  I,  movement  along 


a  hysteresis  curve  3-4,  movement  which  crosses  from  a  hysteresis  to  a  cyclic  curve  1-2,  and  move¬ 
ment  which  crosses  from  one  hysteresis  curve  to  another  hysteresis  curve  6-3.  Using  these  four 
movements  any  complex  response  can  be  built  up  using  Neuber  control. 

Figure  II  shows  Neuber  control  used  to  determine  the  end  point  of  a  local  stress-strain 
movement  corresponding  to  half  a  nominal  load  cycle.  Because  of  the  nature  of  the  cyclic  and 
hysteresis  curves  the  solution  B  must  be  found  by  iteration. 


a 


(3)  Movement  from  hysteresis  to  cyclic  curve  e.g.  1  to  2 

(4)  Movement  from  hysteresis  to  hysteresis  curve  6  to  3 

FIG.  13  CHARACTERISING  LOCAL  STRESS  STRAIN  MOVEMENT 


II 


4.1.3  Numerical  Methods 


Newton’s  method  has  been  used  in  the  implementing  program  since  it  combines  speed  with 
stability  for  functions  like  (2)  and  (5).  The  initial  estimate  of  local  stress  is  the  exact  solution  of 
Neuber  Control  rule  with  the  hysteresis  equation  without  the  elastic  strain  term.  Thus  Ao 
initially  follows  from 


where  TV  1  for  cyclic  curve  (type  of  curve:  7V). 

7V  2  for  hysteresis  curve  (type  of  curse:  TV). 

The  iteration  to  the  solution  from  this  initial  estimate  is  quite  rapid,  even  to  high  accuracy. 
There  will  be  an  error  in  determining  point  B  in  Figure  1 1  regardless  of  the  simulation  techniques, 
so  that  error  control  must  be  a  feature  of  any  Neuber  based  local  stress-strain  simulation. 

In  situations  such  as  depicted  in  Figure  12  where  a  solution  to  an  initial  hysteresis  curve  and 
the  Neuber  control  law  is  B'  the  Newton  iteration  does  not  stop.  Since  B'  is  found  to  exceed 
point  P  then  a  crossover  is  known  to  have  occurred  and  the  iteration  starts  again  using  B'  as 
an  estimate  to  the  solution  which  is  known  to  lie  on  the  previous  hysteresis  curve.  In  complex 
memory  behaviour  the  situation  depicted  in  Figure  12  will  occur  repetitively  with  each  success¬ 
ive  solution  used  as  an  estimate  for  a  further  solution  based  on  a  more  appropriate  curve.  This 
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is  shown  in  Figure  14  and  although  it  may  appear  to  involve  extended  and  complex  calculations 
the  whole  iteration  to  the  correct  solution  is  fairly  direct  within  the  program.  (The  hysteresis 
curves  are  expressed  as  general  functions  that  represent  also  the  cyclic  curve  as  well  as  offset 
situations,  e.g.  the  hysteresis  curve  OB  in  Figure  12  must  be  defined  with  respect  to  the  origin 
A;  the  beginning  of  the  current  movement  in  the  local  stress-strain  stale  )  It  should  be  pointed 
out  here  that  although  there  may  be  several  sets  of  iterations  involved  in  determining  a  solution 
there  is  actually  only  one  inherent  error  step.  This  is  because  intermediate  solutions  such  as 
6"  and  6'  in  Figure  14  are  only  required  to  determine  the  correct  hysteresis  curve.  Only  the 
accuracy  of  the  last  solution  such  as  6.  will  ligurc  in  the  overall  stability  of  the  simulation  during 
the  processing  of  long  histories. 
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FIG.  15  FLOW  CHART  OF  LOCAL  STRESS-STRAIN  SIMULATION 


The  flow  chart  of  the  simulation  by  the  techniques  described,  is  shown  in  Figure  15;  and  its 
implementation  as  a  computer  program  is  discussed  in  the  Appendix. 
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Figure  16  shows  a  portion  of  a  local  stress-strain  history  simulated  in  the  manner  outlined 
from  a  service  load  history  (nominal  strain)  measured  close  to  a  critical  stress  raiser.  Hysteresis, 
memory  and  softening  are  evident,  though  as  previously  mentioned  the  capability  for  the  last 
was  subsequently  considered  unnecessary  and  removed  from  the  model. 


FIG.  16  A  LOCAL  STRESS-STRAIN  HISTORY  SIMULATED  FROM  NOMINAL  STRAIN 
DATA  USING  THE  NEUBER  RULE 


4.2  Miscellaneous  Effects 

In  the  Appendix  it  will  be  seen  that  small  changes  in  the  solution  for  («.  a)  of  (1)  or  (2) 
and  (4)  or  (5)  can  lead  to  quite  different  program  decisions.  However,  in  a  range-pair  sense  the 
differences  arising  thereby  have  only  the  same  order  of  magnitude  as  these  errors,  this  corres¬ 
ponding  to  different  choices  of  turning  points  at  nearly  the  same  stresses. 

In  a  similar  way  the  accumulated  error  in  a  long  sequence  will  most  probably  correspond 
to  an  erroneous  mean  load  for  the  later  range-pairs.  This  accumulated  error  may  also  appear 
on  one  of  the  large  amplitude  range-pairs.  The  reason  for  both  these  effects  is  that  the  conditions 
for  the  smaller  range-pairs  are  likely  to  be  satisfied  over  relatively  few  turning  points  while  the 
larger  range-pairs  may  well  require  the  whole  sequence  for  their  completion.  The  error  in  range 
caused  by  drift  obviously  depends  on  the  number  of  turning  points  between  the  defining  peaks 
and  this  will  probably  be  larger  for  large  range-pairs. 

The  final  solution  of  Equations  (I)  or  (2)  and  (4)  or  (5).  through  the  range-pair  type  of 
algorithm,  depends  uniquely  on  the  immediately  prior  stress-strain  state  and  another  state  further 
in  the  past.  Thus  each  turning  point  is  determined  by  two  others.  Intuitively  one  may  imagine 
that  satisfaction  of  range-pair  conditions  always  corresponds  to  a  closed  hysteresis  loop  in  the 
t-c  plane.  However,  the  mere  dependence  of  each  turning  point  upon  two  previous  ones  is  not 
enough  to  guarantee  this  and  some  examples  without  closure  have  been  found. 
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5.  CALCULATION  OF  FATIGUE  DAMAGE 


Having  translated  a  nominal  load  sequence  into  a  local  stress-strain  sequence  such  as  in 
Figure  16  it  remains  to  resolve  the  damaging  influence  of  the  latter.  The  procedure  used  to 
accomplish  this  consisted  of  two  steps : 

(1)  cycle  counting  of  the  local  stress- strain  history; 

(2)  determination  of  fatigue  damage  for  each  cycle  produced  by  ( 1 )  and  subsequent  sum¬ 
mation. 

The  cycle  counting  method  used  to  process  the  local  stress-strain  histories  was  the  range- 
mean-pair  method-3  which  counts  cycles  that  relate  to  closed  hysteresis  loops.  The  algorithm 
that  constitutes  this  method  is  very  simple,  very  fast  and  counts  in  a  single  pass  which  allows 
it  to  be  run  in  parallel  with  the  simulation  program  (The  simulation  model  actually  incorporates 
range-mean-pair  detection  logic  within  its  own  structure,  though  it  is  easier  to  discuss  the  two 
procedures  separately.  The  logic  is  the  same  as  that  for  selecting  the  proper  hysteresis  curve  in 
solving  for  local  stress  strain  states.) 

Subsequently,  fatigue  damage  for  each  range-mean-pair  cycle  of  local  stress  sir.-on  was 
determined  using  life  data  expressed  as 

A' '  A C.  (6) 

where  N  =  no  cycles  to  failure, 

A  <p  ==  plastic  strain  range  for  local  strain  cycles, 

a,  C  =  material  parameters. 

Equation  (6)  is  commonly  referred  to  as  the  Manson-Coffin  Rule  and  it  represents  the  linear 
relationship  observed  between  the  logarithms  of  the  variables,  i.e. 

log  N  =  or ’(log  C  -  log  A  f,,).  (7) 

Reference  12  contains  cyclic  life  data  for  a  range  of  common  materials  that  allows  estimation 
of  the  constants  Cand  <x  in  (6).  In  general  <x  0-5  0  6  whilst  C  is  a  small  multiple  of  total  strain 
to  static  failure.  The  damage  sustained  per  range-mean-pair  cycle  is  thus 

1/Wj  =  (A  fpJC)'  3,  (8) 

and  for  all  the  cycles  in  a  sequence  1  <  i  <  n  the  total  damage  is 

D  =-  I  "(A  <„,/<: )'  >.  (9) 

I  1 

Hence  by  combining  the  local  stress-strain  simulation  procedure  described  in  the  previous 
section  with  the  range-mean-pair  cycle  counting  technique  and  linear  accumulation  of  Manson- 
Coffin  damage,  a  package  is  obtained  which  can  be  used  for  the  fatigue  analysis  of  any  service 
load  sequence. 

It  should  be  noted  that  though  the  linear  cumulative  damage  law  described  above  proved 
to  be  quite  adequate  in  processing  local  stress-strain  range-mean-pair  cycles,2-5,6-7-8-1318-20  it 
is  thought  that  the  inaccuracies  that  do  exist  in  a  simulation  and  damage  calculation  of  this  sort 
are  due  in  the  main  to  this  linear  rule.8-18  so  that  subsequent  refinements  would  be  best  con¬ 
centrated  in  this  area. 


6.  CONCLUSIONS 

The  advantages  of  utilising  cyclic  stress-strain  concepts  in  the  fatigue  analysis  of  service 
nominal  stress  histories  have  been  well  documented  in  the  literature.  An  algorithm  which  utilises 
these  concepts  in  the  light  of  this  experience  has  been  described  and  its  use  in  the  analysis  of  fatigue 
environments  should  improve  predictions  of  initial  life.  Its  flexibility  and  fast  analysis  of  the 
effects  of  mission  changes,  product  improvements,  load  sequence,  etc.,  should  also  be  of  value. 
It  differs  from  similar  programs  by  finding  numerical  solutions  to  the  original  Neuber  and  stress- 
strain  or  hysteresis  relations  without  approximating  the  model.  The  logical  steps  in  pairing  a 
Neuber  curve  and  a  hysteresis  curve  for  finding  a  stress-strain  state  transpire  to  be  just  the  same 


as  those  used  in  the  stack  inspection  procedure23  for  range-pair  counting.  Since  these  are  merely 
logical  steps  they  cannot  impair  numerical  accuracy.  Because  of  this  and  because  the  model  is 
not  approximated  the  above  method  is  particularly  suitable  for  general,  long,  practical  load 
sequences. 

We  have  already  noted  that  the  application  of  Neuber  analysis  to  the  determination  of 
local  stress-strain  was  such  that  subsequent  fatigue  damage  analysis  was  relevant  to  crack 
initiation  only,  and  that  where  crack  propagation  formed  a  large  portion  of  total  life  an  extended 
analysis  should  be  used  to  incorporate  crack  growth.  The  use  of  a  strain-based  intensity  factor 
solution  for  the  initiation  and  growth  of  short  cracks  from  notches  is  detailed  by  Haddad  et  alP 
an  i  by  combining  this  with  local  stress-strain/damage  analysis  as  described  in  the  present  paper, 
considerable  scope  exists  for  developing  a  general  fatigue  program. 
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FIG.  17  OPERATION  OF  LOCAL  STRESS-STRAIN  STACKS 
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APPENDIX 
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Computer  Simulation  of  Local  Stress-Strain 

Because  of  its  powerful  decision  constructs  the  language  PASCAL  was  used  for  the  program 
of  the  local  stress-strain  simulation  flow-charted  in  Figure  15.  Denoted  LSSF  (Local  Stress- 
Strain  Follower)  the  program  was  divided  into  four  modules: 

(1)  the  function  CURVE  which  determined  the  strains  on  the  hysteresis  (or  cyclic)  and 
Neuber  curves  corresponding  to  a  given  stress  value  as  passed  to  it  by 

(2)  the  function  NP,  containing  a  procedure  based  on  Newton's  method,  which  upon  receipt 
of  an  estimate  to  the  stress  at  the  solution,  where  the  material  response  curve,  either  the 
hysteresis  or  cyclic  curves,  intersects  the  Neuber  curve)  iterates  towards  the  solution 
until  a  specified  relative  error  is  obtained ; 

(3)  the  procedure  MOVE  which  controls  the  output  of  results,  i.e.  the  local  stresses  and 
strains  corresponding  to  input  nominal  stresses  and  strains; 

(4)  the  main  program  which  tracks  the  local  stress-strain  path  incorporating  the  basic 
movements  as  depicted  in  Figure  13. 

The  fundamental  operator  of  the  range-mean-pair  method23 — the  turning  point  stack — was 
used  in  this  program  for  similar  reasons.  A  method  of  storing  local  stress-strain  turning  points 
was  needed  so  that  decisions  about  "cross-overs"  could  be  made.  Figure  9  indicates,  theoretically 
at  least,  that  every  local  stress-strain  turning  point  in  a  history  may  have  to  be  retained  through¬ 
out  a  simulation;  when  the  local  stress-strain  state  moves  in  response  to  the  load  change  6  to  7 
it  jumps  out  of  three  loops  and  information  on  each  of  these  loops  (i.e.  points  6,  5.  4,  3.  2,  1)  is 
required  to  currently  determine  the  point  7.  Thus  stacks  SV[K],  EV[Kj  and  CTYPE[K]  denoting 
respectively  local  stress  stack  of  K  elements,  local  strain  stack  of  K  elements  and  the  type  of 
movement  to  SV[K]  and  EV[K]  (“C"  for  cyclic  curve  and  “H"  for  hysteresis  curve)  were  used 
to  contain  each  local  stress-strain  turning  point  and  the  nature  of  the  movement  involved. 
When  closed  hysteresis  loops  are  detected  the  corresponding  elements  of  the  stacks  are  emptied 
and  the  gaps  closed.  In  this  way  the  stacks  are  continually  loaded  and  emptied  in  response  to 
local  stress-strain  movements.  This  may  be  recognized  as  an  exact  analogy  of  the  operation  of 
the  turning  point  stack  of  Topper  and  Conle13. 

Figure  17  depicts  a  simple  local  stress-strain  movement  using  the  stacks.  The  third  stack 
shows  the  first  approximation  to  the  stress-strain  point  SV[4],  EV[4],  corresponding  to  o'0, 
*V  After  the  third  stack  a  hypothetical  fourth  stack  has  been  added  which  shows  the  values 
oe  and  te,  the  solution  of  correct  equations  derived  from  previous  history  and  the  first  approx¬ 
imations  (see  Figure  18).  In  the  fifth  stack  the  now  established  hysteresis  loop  <n,  o>  has  been 
removed.  In  practice  the  program  proceeds  directly  from  the  third  stack  to  the  fifth. 


The  main  decision  of  the  program  upon  calculation  of  a  solution  point  is  to  determine  if  a 
“cross-over”  occurs  and  if  so  which  type.  This  reduces  to  establishing  whether  the  previous  point 
lies  on  either  a  cyclic  or  hysteresis  curve  and  then  whether  the  present  point  has  crossed  to  the 
outside  of  this  previous  curve.  Figure  IS  illustrates  the  total  simulation  broken  down  into  all 
possible  movements  that  occur  within  the  scope  of  the  concepts  outlined  previously.  It  should  be 
noted  that  although  absolute  stresses  and  strains  are  ultimately  determined  the  solution  of  a  given 
material  response  using  Neuber,  is  determined  with  respect  to  the  immediately  preceding  local 
stress-strain  state,  i.e.  the  local  stress-strain  co-ordinates  implied  by  the  Neuber  equation  as 
depicted  in  Figure  1 1  are  used.  Decisions  about  equality  of  real  numbers  must  be  made  with 
regard  to  a  relative  error  criterion,  e.g.  two  stresses  or  strains  may  be  considered  equal  if  they 
differ  by  less  than  0-005°o. 

0  Consider  the  local  stress-strain  response  that 

i  1  may  result  from  the  Kr/i  nominal  load  application 

/  AS  (the  prior  local  history  has  been  determined  and 

/  is  contained  in  the  stacks  SV.  EV,  CTYPE). 


(I)  An  estimate  of  the  res¬ 
ponse  is  made  using  the 
material  hysteresis  curve. 
The  stacks  are  altered 
accordingly : 
i.e.  SV[KJ  =  o' 
EV[K]  =  t' 
CTYPE[K]  -  “H” 


Ao  Ae  =  (KtASl2/E 


(2)  This  response  is  checked  against  prior  material  behaviour  using  CTYPE. 

For  CTYPE  [K  —  1]  =  “C",i.e.  previous  point  lies  on  cyclic  curve,  three  movement  types 
are  possible:  >  \  \  a 
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|SV[K]|  -  |SV[K  -  1]| 
SV[K]  =  o' 

EV(K]  =  t' 
CTYPEfK]  -  “H" 


;sv[K];  -  sv[K  -  i]  iSV[K]  sv(k  -  iji 

SV[K]  -  o'  SV[K]  =  a” 

EV[K]  ■  «'  EV[K]  - 

CTYPE[K]  "C”  CTYPE(K]  *C" 


(Crossover  from  hysteresis  curve  must  occur.  Correct  point  o'\  t"  determined  using  cyclic  curve.) 


For  CTYPEfK  —  1]  ^  “C”,  that  is,  the  previous  point  lies  on  a  hysteresis  curve,  the 
resulting  movement  will  depend  on  SV[K  —  2]  and  CTYPEfK  -  2].  For  example  let 

TEST  -  ;SV[K ]  SVfK  -  1]  . 


Then  for  TEST  =  0  the  movement  is  such  that 


SV[K]  =  a', 
EVfK]  =  «\ 
CTYPE[K]  -  “H“. 


SVfK  —  t] 
EVfK  -  1] 
CTYPEfK-  11 


SV1K-2J 


Otherwise  for  TEST  >  0  two  movements  are  possible: 

(i)  CTYPEfK  -  2]  =  "C” 

Stress-strain  state 
moves  from  hysteresis 
curve  to  cyclic  curve 
as  shown: 

SV[K]  =  o', 

EV[K]  = 

CTYPEfK]  =  “C”. 


SVIK-1] 
EV[K-  11 
CTYPEfK  -  1]  ='H' 


SVIK-21 
EV[K  —  21 
CTYPEfK -2]  =4C' 


(ii)  CTYPEfK  -  2]  *  “C”. 

Stress-strain  state  moves  from  hys¬ 
teresis  curve  to  hysteresis  curve. 

For  complex  memory  situations  in¬ 
volving  a  movement  out  of  many 
loops  the  local  state  may  eventually 
jump  out  of  cyclic  curve,  i.e.  move¬ 
ments  according  to  (i)  above  only 
occur  at  the  end  of  this  movement. 

Note  that  contents  of  stack  are 
overwritten  as  in  Figure  17  when 
loops  are  closed  by  a  movement  in 
state. 


FIG.  18  PROGRAM  DECISIONS  AND  RESULTING  LOCAL  STRESS-STRAIN  STATE 
MOVEMENTS  (contd.) 
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